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Abstract. We report new follow-up spectroscopy of i'-dropout galaxies with an NB921-band depression found 
in the Subaru Deep Field. The NB921-depressed i'-dropout selection method is expected to select galaxies with 
large equivalent width Lya emission over a wide redshift range, 6.0 z S3 6.5. Two of four observed targets show 
a strong emission line with a clear asymmetric profile, identified as Lya emitters at z = 6.11 and 6.00. Their rest- 
frame equivalent widths are 153A and 114A, which are lower limits on the intrinsic equivalent widths. Through 
our spectroscopic observations (including previous ones) of NB921-depressed i'-dropout galaxies, we identified 
5 galaxies in total with a rest-frame equivalent width larger than IOOA at 6.0 ^ z 6.5 out of 8 photometric 
candidates, which suggests that the NB921-depressed i'-dropout selection method is possibly an efficient way to 
search for Lya emitters with a large Lya equivalent width, in a wider redshift range than usual narrow-band 
excess techniques. By combining these findings with our previous observational results, we infer that the fraction 
of broad-band selected galaxies having a rest-frame equivalent width larger than 100A is significantly higher at 
z ~ 6 (the cosmic age of ~1 Gyr) than that at z ~ 3 (~2 Gyr), being consistent with the idea that the typical 
stellar population of galaxies is significantly younger at z ~ 6 than that at z ~ 3. The NB921-depressed i'-dropout 
galaxies may be interesting candidates for hosts of massive, zero- met allicity Population III stars. 

Key words. early universe - galaxies: evolution - galaxies: formation - galaxies: individual (SDF 
J132345. 6+271701, SDF J132519.4+271829) - galaxies: starburst 



1. Introduction 

Observational searches for high-z galaxies have achieved 
important progress in this decade. Thanks to large tele- 
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Table 1. Photometric properties of NB921-depressed z'-dropout galaxies in the Subaru Deep Field 
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SDF J132521. 6+274229 
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24.73 
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a AB magnitude with a 2" aperture photometry. Lower limits are 2a limiting magnitudes. 
b References. — 1: This work, 2: Nagao et al. (2005a), 3: Nagao et al. (2004). 
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Fig. 1. Color-magnitude diagram of z 1 — NB921 vs. z' 
{left) and z' — NB921 vs. i' — z' color-color diagram 
(right). The 48 z'-dropout objects detected in the SDF 
are plotted. The NB921-depressed z'-dropout objects are 
shown by open squares, and the other z'-dropout objects 
are shown by filled circles. Error bars denote lc uncer- 
tainty in the observed color or magnitude. In the left 
panel, 1, 2, and 3 <r uncertainties of the z' — NB921 color 
are shown by the dotted, dashed, and solid lines, respec- 
tively. 



scopes such as Subaru, VLTs and Keck, a few dozens of 
galaxies at 6 < z < 7 have been spectroscopically con- 
firmed so far (Hu et al. 2002; Kodaira et al. 2003; Cuby 
et al. 2003; Rhoads et al. 2004; Kurk et al. 2004; Nagao 
et al. 2004, 2005a; Taniguchi et al. 2005; Stern et al. 2005; 
Kashikawaet al. 2006). Statistical properties of these high- 
z galaxies, such as their luminosity function and their 
correlation function provide information on galaxy evolu- 
tion, the cosmic star-formation rate, and the re-ionization 
history of the universe. In addition to statistical proper- 
ties, properties of individual galaxies at high redshift also 
shed light on galaxy evolution. In particular, the equiva- 
lent width (EW) of the Lya emission provides clues on 



the stellar population of high-z galaxies, for which it is 
generally difficult to investigate stellar spectral features 
due to the limited observational sensitivity of currently 
available instruments. Malhotra & Rhoads (2002) pre- 
sented their theoretical estimates that large Lya equiv- 
alent widths [EWq ^ 150 A for a Salpeter initial-mass 
function (IMF) and EW £ 240A for a flatter IMF] sug- 
gest a very young (< 10 7 years) stellar population (e.g., 
Tumlinson et al. 2003). More interestingly, it has been the- 
oretically predicted that galaxies hosting zero-metallicity 
stars (or Population III stars; hereafter PopIII) should 
show huge Lya equivalent widths which could reach up 
to EW a ~ 1000A (e.g., Schaerer 2002, 2003; Tumlinson 
et al. 2003; Scannapieco et al. 2003). Therefore, the fre- 
quency distribution of the Lya equivalent width in high-z 
galaxies and its evolution with redshift are important in 
understanding the early stages of galaxy evolution (e.g., 
Shimasaku et al. 2006; Ando et al. 2006). 

However, the observational study of the Lya equivalent 
width in galaxies at high redshift is not straightforward. 
Although Malhotra & Rhoads (2002) reported that more 
than half of their sample of 150 Lya emitters (LAEs) have 
Lya equivalent width larger than 240A, their analysis is 
based only on photometric data. Especially for LAEs at 
z ~ 6.5 (which corresponds to the window between airglow 
emission at A ~ 9200 A), the z'-band magnitude is con- 
taminated by the Lya emission, and thus an unreliable 
measure of continuum flux density. Even for the sample 
with spectroscopic data, the continuum emission of LAEs 
selected by using narrow-band magnitude is generally too 
faint to be detected in their spectra. 

We are exploiting a new selection method to iden- 
tify LAEs with a large Lya equivalent width, that is, 
z' dropout with a "depression" in the narrow-band filter 
NB921 (see §2) (Nagao et al. 2004, 2005a). If galaxies are 
at 6.0 < z < 6.5, their redshiftcd Lya emission is ex- 
pected at 8500A ^ A b s ^ 9100A. In this case, the narrow- 
band magnitude at A ~ 9200A has no contamination from 
Lya emission and thus traces only the continuum, while 
the z'-band magnitude is enhanced by the contribution of 
the strong Lya emission and thus the narrow-band mag- 
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nitude is depressed with respect to the z'-band magni- 
tude. Therefore this selection method finds galaxies with 
a large Lya equivalent width in a wide redshift range, 
6.0 < z < 6.5. Note that this selection method does not 
select Galactic late- type stars, unlike usual z'-dropout se- 
lection (see Nagao et al. 2005a for more details). 

In this paper, we report new spectroscopy of a sample 
of NB921-depressed z'-dropout galaxies. Throughout this 
paper, we adopt a cosmology with (Qtot, ^m, ^a) = (1-0, 
0.3, 0.7) and H = 70 km s" 1 Mpc" 1 . 

2. Target selection and observation 

We selected our z'-dropout galaxy sample from the pub- 
lic catalog of the Subaru Deep Field (SDF) imaging sur- 
vey (Kashikawa et al. 2004), which contains broad-band 
(B, V, Rq, z', and z') and narrow-band photometric data 
[NBS1Q and A^_B921; the central wavelengths and the 
half-widths of the transmittance are (8150A, 120A) and 
(9196A, 132A), respectively]. The adopted criteria to se- 
lect the z'-dropout galaxy sample are: 

— z' < 26.1 (i.e., above 5a background error), 

— i' -z' > 1.5, 

— B > 28.5 (below 3a background error), and 

— Rc > 27.8 (below 3a background error). 

Among the 48 selected z'-dropout galaxies in SDF, we 
identified as a "NB921-depressed z'-dropout galaxy" an 
object whose z' — NB921 is less than — 2a of the sky noise. 
Here "sky noise" refers to the error on the z' — NB921 
color due to sky noise in the z' and NB921 filters. The 
z' — NB921 versus z' color-magnitude diagram and the 
z' — NB921 versus i' — z' color-color diagram are shown 
in Figure 1. Eight z'-dropout galaxies satisfy this crite- 
rion and thus are classified as NB921-depressed z'-dropout 
galaxies. Note that the significance level of the NB921 de- 
pression of one of the eight NB921-depressed z'-dropout 
galaxies (SDF J132440.6+273607; Nagao et al. 2004) is 
slightly less than 2a. However we included this object in 
our NB921-depressed z'-dropout galaxy sample. The ob- 
ject names and the photometric properties of this sample 
are summarized in Table 1, and the thumbnails of these 
objects are shown in Figure 2. 

We have already obtained optical spectra 
of 3 NB921-depressed z'-dropout galaxies; SDF 
J132440.6+273607 (#4; z = 6.33; Nagao et al. 2004), 
SDF J132426.5.0+271600 and SDF J132442. 5+272423 
(#3 and #5; z = 6.03 and 6.04; Nagao et al. 2005a). 
Among the remaining 5 NB921-depressed z'-dropout 
galaxies, we observed 4 objects on 26 April 2006 (UT). 
This spectroscopy was carried out with the Faint Object 
Camera And Spectrograph (FOCAS; Kashikawa et al. 
2002) on the Subaru telescope (lye et al. 2004), in its 
multi-object slit mode. The four objects were chosen to 
optimize the slit mask design (by also accomodationg 
several other objects observed within the broader con- 
text of the SDF project). The 175 lines mm -1 echelle 
grating and the SDSS z' filter were used. The resulting 
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Fig. 2. Thumbnail images of NB921-depressed z'-dropout 
galaxies. The square regions around each object in the z', 
z', and NB921 images are shown from left to right. The 
IDs given at the left side of the panels correspond to those 
in Table 1. The size of panels and the radius of circles are 
16 arcsec and 8 arcsec, respectively. 



wavelength coverage was ~ 8300 — 10000A, with a 
dispersion of ~ 0.95 A pixel -1 . The adopted slit width 
was 0.83 arcsec, giving a spectral resolution of R ~ 1500 
or AA ~ 6A at ~ 9000A as measured from the widths 
of atmospheric OH emission lines. The spatial sampling 
was 0.31 arcsec per resolution element, as we adopted 
3 pixel on-chip binning. The seeing was variable during 
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Table 2. Spectroscopic properties of NB921-depressed i'-dropout galaxies in the Subaru Deep Field 



No. 


Redshift 


F(Lya) a 


FWHM b 


EW (Lya) c 


£W / min (Lya) d 


Ref. c 






(10~ 17 ergs/s/cm 2 ) 


(km s- 1 ) 


(A) 


(A) 
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6.11 


4.5±0.1 


330 


153 


72 


1 


3 


6.03 


3.6±0.3 


410 


94 


83 


2 


4 


6.33 


4.0±0.1 
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130 


133 


3 
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6.04 


4.5±0.3 


220 


236 


146 


2 


6 


6.00 


3.4±0.2 


220 


114 


90 


1 



a No correction for the slit loss. 

b Corrected for the instrumental broadening (FWHM = 6.2A). 

c Not corrected for the absorption effects due to the Lya-forest, and thus lower limits on the intrinsic equivalent width. 

d Minimum _EWo(Lya) that can be selected by the NB921-depressed i'-dropout method for galaxies with the redshift and the 

z'-band magnitude of each NB921-depressed i'-dropout galaxy, estimated by assuming a flat UV slope (see §§4.1). 
e References. — 1: This work, 2: Nagao et al. (2005a), 3: Nagao et al. (2004). 
f Observed line width is comparable to the instrumental broadening width. 



the observation (0.5 - 1.0 arcsec). The total integration 
times were 7200 sec for SDF J132345. 6+271701 (#1) 
and SDF J132422.0+271742 (#2), and 9000 sec for SDF 
J132519.4+271829 (#6) and SDF J132526. 1+271902 
(#8). Note that the position angle of the slits was set 
to the North-South direction for all the target objects. 
Since no neighboring objects are around the target 
objects (Figure 2), the detected signal is believed to 
be from the targets. We also obtained spectra of the 
spectrophotometric standard star Feige 34 (Oke 1990) for 
flux calibration. The obtained data were reduced in the 
standard manner, by using IRAF. 

3. Results 

Among the four observed objects, SDF J132345. 6+271701 
(#1) and SDF J132519.4+271829 (#6) show a strong 
emission line in their spectra, whose peak wavelengths 
are at 8634A and 8512A, respectively. In Figures 3 and 

4, the sky-subtracted position-velocity spectrogram, the 
extracted one-dimensional spectrum, and the typical sky 
spectrum are shown for these two objects, respectively. 
The aperture size used for the extraction of the one- 
dimensional spectrum is 5 binned pixels (~1.6 arcsec). 
No continuum emission is detected for either object. Both 
emission lines show a clear asymmetric profile, i.e., a sharp 
decline on the blue side and a prominent tail on the red 
side, which suggests that the observed emission lines are 
Lya. Note that if the detected emission line were H/3, [O 
iii]A5007 or Ha, other rest-frame optical emission line(s) 
should be seen in the observed wavelength range. If the de- 
tected emission line was [O n]A3727, it should be resolved 
as a doublet emission, since the expected wavelength sep- 
aration of the redshifted [O n]AA3726,3729 is - 6 A (or 
~6 pixels), which corresponds to ~ 210 km s _1 . Here 
we assume that the velocity width of the [O n] lines is 
not very broad; otherwise we could not resolve the [O II] 
doublet. However our assumption seems valid for normal 
star-forming galaxies, because the [O n] doublet of some 



star-forming galaxies in the SDF is indeed resolved in our 
previous spectroscopic follow-up observations (Shimasaku 
et al. 2006; Ly et al. 2007). The fact that we see only one 
strong emission line with an asymmetric profile in each 
spectrum strongly suggests that we have detected Lya at 
2~ 6. 

To quantify the asymmetry of the detected emis- 
sion lines, we calculated two independent parameters; 
/red //blue and S w . The former is the ratio between / rcc i 
and /blue, where / ro d is the flux at wavelengths longer 
than the emission-line peak, while /blue is that at shorter 
wavelengths (Taniguchi et al. 2005; Nagao et al. 2005a; 
see also Haiman 2002 for a theoretical discussion on the 
parameter / re d//biue)- The measured ratios are 2.08+0.12 
and 2.47+0.30 for SDF J132345.6+271701 (#1) and SDF 
J132519.4+271829 (#6), respectively. Note that the lat- 
ter value may be overestimated, since the OH airglow 
emission at the blue side of the Lya emission may be 
over-subtracted (Fig. 4). The / rc d//biue ratio of both line- 
detected objects is significantly larger than unity, and also 
larger than that of most LAEs reported by Taniguchi et 
al. (2005). This result is consistent with the interpreta- 
tion that the detected emission lines are Lya. The lat- 
ter parameter (S w ) is the weighted skewness (Shimasaku 
et al. 2006; Kashikawa et al. 2006), which is larger for 
objects with higher asymmetries and/or larger emission- 
line widths. The calculated values are 4.00+0.22 and 
8.72+0.28 for SDF J132345. 6+271701 (#1) and SDF 
J132519.4+271829 (#6), respectively. Since the weighted 
skewness of emission lines arising from possible low- 
redshift interlopers generally do not exceed 3 (Kashikawa 
et al. 2006; see also Shimasaku et al. 2006), the derived val- 
ues of the weighted skewness are also consistent with the 
interpretation that the detected emission lines are Lya. 
Taking the derived values of / re d//biue and S w for the two 
line-detected objects into account, we conclude that the 
detected emission lines are Lya and accordingly that the 
redshifts of the two objects are 6.11 and 6.00, respectively. 
The observed Lya fluxes are (4.5+0.1) x 10 -17 ergs s -1 
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Fig. 3. Sky-subtracted optical position-velocity spectro- 
gram (top) and one-dimensional spectrum (middle) of 
SDF J132345. 6+271701 (#1) obtained with FOCAS on 
Subaru. The typical spectrum of the sky emission is also 
shown in the bottom panel. 



cm~ 2 and (3.4±0.2) x 10 -17 ergs s _1 cm~ 2 , without any 
correction for slit loss. The Lya luminosities are thus cal- 
culated to be 1.9 x 10 43 ergs s -1 and 1.3 x 10 43 ergs s _1 , re- 
spectively. The measured emission-line widths in FWHM 
are 11.3l}'gA and 8.6i~g'yA. These values correspond to 

9.6l2:o A and Q - 2 -H^, or velocity widths of 330t™ km 
s~ 4 and 220^220 km s _1 , respectively, after corrected for 
the instrumental broadening effect. Here the absorption 
effects are not taken into account. 

Adopting these redshifts, the NvA1240 emission line 
would be expected at 8816A and 8680A for SDF 
J132345.6+271701 (#1) and SDF J132519.4+271829 
(#6) respectively, if these two objects were active galactic 
nuclei (AGNs). However, no emission-line features are seen 
at the corresponding wavelengths (see Figures 3 and 4). 
The 3a upper limits for the N v emission are 5.7 x 1CP 18 



ergs s 1 cm 2 and 6.9 x 10 



18 ergs s 1 cm 2 for these two 



objects respectively, assuming the same velocity widths as 
Lya. However, we cannot completely rule out the possibil- 
ity that one or both of these two objects are AGNs, since 
some high-z narrow-line radio galaxies show very weak N V 
compared with Lya (e.g., De Breuck et al. 2000; Nagao et 
al. 2006; see also Malkan et al. 1996). 

As mentioned in §1, it is sometimes very difficult to 
estimate the Lya equivalent width because measuring 
the continuum in the low-S/N spectra is very difficult. 
However, we can determine the flux density of the contin- 
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Fig. 4. Same as Figure 3 but for SDF J132519.4+271829 
(#6). 



uum emission for those two galaxies from the NB921-band 
magnitude, since the NB921-band flux does not contain 
Lya flux, unlike the case of LAEs selected by the NB921 
excess (e.g., Kodaira et al. 2003; Taniguchi et al. 2005). 
The NB921-band magnitudes of the two galaxies are 26.37 
mag and 26.38 mag, which correspond to the flux densi- 
ties of fx = (4.1 ± 1.1) x 10~ 20 ergs s _1 cm~ 2 A -1 and 
-20 ___.„ „-i 



ergs s 

(4.2±1.2) x 10~ 20 ergs s _1 cm~ 2 A at the wavelengths of 
the Lya peak, respectively. Then the rest-frame equivalent 
widths of the Lya emission are calculated to be 153 ±42 A 
and 114±32Afor SDF J132345. 6+271701 (#1) and SDF 
J132519.4+271829 (#6) respectively, assuming a flat UV 
continuum. Note that if the detected emission lines were 
[O 11], the rest-frame equivalent widths would be 470A 
and 349A, respectively. These large values are quite rare 
for [O 11] emitters (e.g., Ajiki et al. 2006), which also sup- 
ports our interpretation that the detected emission lines 
are Lya. 

The obtained spectroscopic properties are summarized 
in Table 2. In the same table, properties of our previous 
spectra (Nagao et al. 2004, 2005a) are also given for the 
reader's convenience. Note that there are no spectral fea- 
tures in the obtained spectra of SDF J132422.0+271742 
(#2) and SDF J132526. 1+271902 (#8). We will discuss 
these non-detections briefly in the next section. In this 
table, we also give the minimum £Wo(Lya) that can be 
detected by the NB921-depressed i'-dropout method for 
galaxies with the redshift and the z'-band magnitude of 
each NB921-depressed i'-dropout galaxy. We will discuss 
this quantity in §§4.1. 
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4. Discussion 

4.1. The success rate of the NB921-depressed 
i-dropout selection method 

In this observing run we observed four NB921-depressed 
z'-dropout galaxies, and found strong emission lines in 
two of them. Could the two remaining emission-line un- 
detected objects also be strong LAEs at 6.0 < z < 6.5, 
as expected by the NB921-dcprcssed z'-dropout selection? 
Since high-z Lya emission shows a clear asymmetric pro- 
file and especially a prominent tail toward the red side 
of the emission-line peak, the Lya emission of LAEs is 
expected to be resolved with the current wavelength reso- 
lution R ~ 1500. This means that the Lya emission should 
be found even when the redshiftcd Lya line falls on an iso- 
lated (i.e., unresolved) OH airglow emission line. However, 
the Lya detection would be difficult when the Lya emis- 
sion line falls on blended OH lines, such as the lines at 
8760A <, A ^ 8780A seen at the bottom panel of Figure 3 
(which corresponds to a redshift range 6.21 < z ^ 6.22). 
SDF J132422. 0+271742 (#2) has z' and NB921 magni- 
tudes which are the faintest among the NB921-depressed 
z'-dropout galaxy sample. Thus the non-detection of the 
emission line might be simply due to an insufficient inte- 
gration time. On the contrary, the z' and iVi3921 mag- 
nitudes of SDF J132526. 1+271902 (#8) are the brightest 
among the NB921-depressed z'-dropouts. Since this ob- 
ject is so bright, it is selected as an NB921-depressed ob- 
ject in spite of its small amount of the NB921 depression, 
z' — NB921 = —0.25. Therefore the expected Lya equiva- 
lent width of this object is so small (which is estimated to 
be ~ 50A in the rest frame by assuming z = 6.0 and a flat 
UV slope) that the detection of Lya could be more diffi- 
cult than the other NB921-depressed z'-dropout galaxies. 
Taking all of these considerations into account, we can- 
not reject the possibility that the two Lya non-detected 
objects are also LAEs at 6.0 < z < 6.5. Deeper spec- 
troscopic observations are necessary to investigate these 
objects further. 

Therefore, the "success rate" of the NB921-depressed 
z'-dropout selection method (i.e., the probability that this 
selection identify LAEs at 6.0 < z < 6.5 correctly) is 
between 5/8 (63%) and 8/8 (100%). The rest-frame Lya 
equivalent width of all of the five Lya-detected objects 
is quite large, at least ~ 100A, as also expected by the 
selection method (see Table 2). It should be mentioned 
here that, as discussed in Nagao et al. (2005a), the NB921 
depression is not expected in Galactic late-type stars and 
thus this selection method is not contaminated by stars, 
unlike a simple z'-dropout sample (see, e.g., Stanway et al. 
2004). 

Here we discuss on the limiting EW(Lyot) of the 
NB921-depressed z'-dropout selection method. This selec- 
tion technique does not select only LAEs with a large 
EWo(Lya) in principle, because the selection criterion on 
the NB921 depression is based on the sky error as de- 
scribed in §2. This is different from usual LAE surveys 
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Fig. 5. Minimum EWo(Lya) that can be selected by the 
NB921-depressed z'-dropout method [£W min (Lya)] as a 
function of the z'-band magnitude, estimated by assuming 
z = 6.0 and a flat UV slope. 



based on the narrow-band excession, for which a certain 
excess magnitude is generally used instead of the sky noise. 
Therefore the NB921-depressed z'-dropout method would 
select LAEs with a relatively small EWq (Lya) if the imag- 
ing data were enough deep and accordingly the sky noise 
was very small. In other words, the minimum EWo(Lya) 
that can be selected by the NB921-depressed z'-dropout 
method [hereafter EW™ 111 (Lya)] depends on the magni- 
tude of target galaxies, since the z' — NB921 color of 
brighter objects is less affected by the sky error with re- 
spect to that of fainter objects. To see this effect quanti- 
tatively, we show the dependence of EW™ m (Lya) on the 
z' magnitude for the case of the SDF dataset in Figure 
5, where z — 6.0 and a flat UV slope are assumed. It is 
demonstrated that only LAEs with EW (Lya) Z 100 A 
are selected among faint galaxies with z' <; 25.5, while 
LAEs with E Wo (Lya) <~ a tens of 10A could be also se- 
lected among relatively bright galaxies with z' ^ 25.0. 
Note that the latter case corresponds to the case of SDF 
J132526. 1+271902 (#8). To compare EW miTl (Lya) and 
the detected EWo(Lya) for each NB921-depressed z'- 
dropout galaxy, we also give EW™ m (Lya) for each set of 
the redshift and the z'-band magnitude of 5 spectroscop- 
ically identified NB921-depressed z'-dropout galaxies in 
Table 2. In order to carry out systematic surveys for strong 
LAEs utilizing the NB921-depressed z'-dropout method, it 
would be more appropriate to adopt a certain z' — NB921 
depression magnitude as a selection criterion rather than 
using the sky error. The reason for adopting the selection 
criterion based on the sky noise is that one of our motiva- 
tions in the photometric selection is to construct a sample 
of the target objects for our spectroscopic observations to 
find candidates of PopIII-hosting galaxies (see Nagao et 
al. 2005b). Due to the limited observing time, we had to 
focus on convincing candidates, i.e., objects with a signif- 
icant NB921 depression with respect to the sky error. 
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4.2. Implication for the stellar population 

The rest-frame Lya equivalent width of the NB921- 
depressed i'-dropout galaxies with a Lya detection (given 
in Table 2) ranges from 94A to 236A. Although these val- 
ues are less than the "critical value" EW (Lya) = 240A 
above which it cannot be explained by normal stellar pop- 
ulations (e.g., Malhotra & Rhoads 2002), it should be 
noted that the values we obtained are not corrected for 
absorption effects. The amount of the Lya absorption de- 
pends on various parameters including the neutral frac- 
tion of the inter-galactic matter (IGM) and the kinematic 
status of neutral hydrogen within LAEs themselves (e.g., 
Ncufeld 1990; Haiman 2002; Mas-Hesse et al. 2003; Ahn 
2004). Here it should be kept in mind that the IGM is not 
perfectly re-ionized at z ~ 6, which is suggested by recent 
observations of high-z quasars (Fan et al. 2006) and LAEs 
(Kashikawa et al. 2006). It has been extensively argued 
whether Lya photons from galaxies at an earlier epoch 
than the completion of the cosmic rc-ionization arc ob- 
servable or not. Some calculations suggest that the Lya 
photons from such galaxies are strongly suppressed (e.g., 
Miralda-Escude 1998; Miralda-Escude & Rees 1998; Loeb 
& Rybicki 1999) and other theoretical works infer that 
a large fraction of the Lya photons can transmit due 
to cosmological H n regions (e.g., Madau & Rees 2000; 
Cen & Haiman 2000). Haiman (2002) investigated vari- 
ous parameter dependences of the Lya transmission frac- 
tion for high-z LAEs by taking both the neutral hydrogen 
within LAEs themselves and the neutral IGM into ac- 
count. They showed that the ratio of the transmitted Lya 
flux to the total (intrinsic) Lya flux [Fo(Lya) / F ioi (Lya)] 
is estimated to be ~10% for typical LAEs with a star- 
formation rate (SFR) - 10 M yr" 1 and -50% for LAEs 
with SFR ~ 100M Q yr" 1 . Following this result, we as- 
sume F (Lya)/F to t (Lya) ~0.5 to correct the absorption 
effect on the observed Lya flux rather conservatively, in 
the sense that the actual Fo (Lya) /F ioi (Lya) of our spec- 
troscopic sample would be smaller than 0.5 since the SFR 
is less than 100 M Q yr" 1 (Nagao et al. 2004, 2005a). See 
also, e.g., Santos (2002) and Dijkstra et al. (2006) for the 
justification of our adopted value of Fo (Lya) /F tot (Lya) . 
Interestingly, by adopting this assumption, the intrin- 
sic £Wo(Lya) for three out of the five spectroscopically 
identified NB921-dcprcsscd i'-dropout galaxies is then ex- 
pected to exceed the critical value £Wo(Lya) = 240A. 
This result is consistent with the idea that galaxies at 
z <; 6 contains young stellar populations, whose age may 
be younger than — 10 7 years as discussed by Malhotra & 
Rhoads (2002); see also Shimasaku et al. (2006). 

The huge intrinsic Lya equivalent width would indi- 
cate that the NB921-depresscd z'-dropout galaxies may 
contain a significant number of PopIII stars. Therefore 
our NB921-depressed i'-dropout galaxy sample offers fas- 
cinating targets for observational searches for PopIII 
stars, which will be an important topic in the forth- 
coming decade (see, e.g., Jimenez & Haiman 2006). We 
have already searched for an observational signature in 



one of the NB921-depressed i'-dropout galaxies, SDF 
J132440.6+273607 (#4; Nagao et al. 2005b), through He 
nA1640 emission (see, e.g., Tumlinson & Shull 2000; Oh et 
al. 2001; Schaercr 2002, 2003). Although the He nA1640 
emission was not detected in this object, similar obser- 
vations will provide important constraints on theoretical 
PopIII models (e.g., Scannapieco et al. 2003; Nagao et al. 
2005b). 

4.3. Evolution of the Ly alpha equivalent width 
distribution 

It has been known that narrow-band selected galaxies 
(e.g., emission- line galaxies) at high rcdshift tend to have 
high Lya equivalent widths, which sometimes exceeds 
the critical value of EW (Lya) = 240A (e.g., Malhotra 
& Rhoads 2002; Shimasaku et al. 2006). However, some 
broad-band selected high-z galaxies such as Lyman-break 
galaxies (LBGs) rarely have high Lya equivalent widths. 
For instance, at z ~ 3 the fraction of LBGs showing 
EW (Lya) > 100A is wl% of -1000 galaxies (Shapley 
et al. 2003). Since our sample is also broad-band selected, 
it is interesting to compare the fraction of galaxies with 
a large equivalent width to investigate whether the stellar 
population of galaxies evolves as a function of redshift. At 
z — 6, our i'-dropout galaxy sample (48 objects) contains 
at least 5 galaxies with £W (Lya) <; 100A [here we re- 
tain SDF J132426. 5+271600 (#3) as a large EW object 
although its EW is slightly below the criterion] . However, 
z'-band flux of the NB921-depressed i'-dropout galaxies is 
strongly boosted by the Lya contamination. After correc- 
tion for this effect, they may fail to satisfy the i'-dropout 
criterion, i' — z' > 1.5. Indeed if adopting the NB921 mag- 
nitude instead of z' as a continuum magnitude at the long 
side of the Lyman break, SDF J132442. 5+272423 (#5) 
should be removed from the SDF i'-dropout galaxy sam- 
ple (i' — NB921 = 0.98). Therefore we estimate the frac- 
tion of i'-dropout galaxies having £Wo(Lya) <; 100A is 
4/48 w 8%. Note that this estimated fraction is a con- 
servative lower limit, because (1) some Lya non-detected 
NB921-depressed i'-dropout galaxies may also have such 
a high equivalent width as discussed above, and (2) some 
of the 48 i'-dropout objects may be emission-line galaxies 
at lower redshift or Galactic late-type stars, since most of 
the SDF i'-dropout objects have not been confirmed by 
spectroscopic follow-up observations (see, e.g., Stanway et 
al. 2004). 

The derived lower limit of the fraction of i'-dropout 
objects having EWo(Lya) il 100A is significantly higher 
than the value at z — 3. However, it should be noted that 
the limiting magnitude of these two samples are different; 
that is, while the spectroscopic survey of LBGs at z — 3 
by Shapley et al. (2003) reaches down to Rab ~ 25.5, our 
SDF i'-dropout sample consists of galaxies with z' < 26.1. 
These limiting magnitudes correspond to the absolute UV 
magnitudes of Mi 50 o = -20.0 and -20.6 respectively, and 
thus the galaxies in the z — 6 sample are intrinsically 
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brighter with respect to those in the z ~ 3 sample, sys- 
tematically. Recently, Ando et al. (2006) reported that 
the LBGs with a high Lya equivalent width arc rarer in 
brighter samples than in fainter samples, at z > 5. Shapley 
et al. (2003) also reported based on their huge number of 
a LBG spectroscopic sample that the broad-band magni- 
tude of LBGs with a stronger Lya emission tends to be 
fainter than that of LBGs with a weaker Lya emission (or 
with a Lya absorption instead of emission; see Table 3 of 
Shapley et al. 2003). Therefore, taking the difference in 
the limiting magnitude between the samples at z ~ 6 and 
at z ~ 3 and the dependence of the equivalent width on 
the luminosity into account, the difference in the fraction 
of galaxies with EWo(Lya) Z 100A between at z ~ 6 and 
at z ~ 3 should be even more significant. Note that the 
significance may be intrinsically much more if the IGM 
effect against the Lya transmission is stronger at z ~ 6 
than at z ~ 3. This result is consistent with a recent study 
by Shimasaku et al. (2006) that the fraction of LBGs hav- 
ing EWo(Lya) Z 100A significantly evolves from z ~ 3 
to z ~ 6. All of these results are consistent with the idea 
that the typical stellar population of galaxies is signifi- 
cantly younger at z ~ 6 than that at z ~ 3. This may 
also be consistent with recent findings that the slope of 
the rest-frame UV continuum of some galaxies is bluer at 
z ~ 6 than the typical UV slope at z ~ 3 (e.g., Stanway 
et al. 2004, 2005; Bouwens et al. 2005; Yan et al. 2005). 

5. Summary 

In earlier papers we identified 8 NB92I-depressed %'- 
dropout galaxies, which are expected to be strong LAEs 
at 6.0 z ^ 6.5, through the narrow-band and broad- 
band photometric data of the SDF (Nagao et al. 2004, 
2005a). In addition to 3 previously spectroscopically con- 
firmed ones, we found that other two NB92I-depressed i'- 
dropout galaxies are also LAEs with EWo(Lya) > 100A at 
z — 6.11 and 6.00, by new optical spectroscopy. This result 
combined with our previous spectroscopic runs means that 
at least 5 objects among 8 NB921-depressed «'-dropout 
galaxies are indeed LAEs having E Wq (Lya) is 100A at 
6.0 ^ z ^ 6.5; these results suggest that the NB921- 
depressed i'-dropout selection method is an efficient tech- 
nique to identify strong LAEs in a wide redshift range, 
6.0 < z < 6.5. 

The obtained results also suggest that more than 8% 
of the «'-dropout galaxies in the SDF have a large Lya 
equivalent width of W (Lya) Z 100A. This is in contrast 
with LBGs at z ~ 3, where such strong LAEs are much 
rarer (~1%). This also implies a strong redshift evolution 
in the Lya equivalent width distribution from z ~ 6 to 
z ~ 3, consistently with a stellar population of broad- 
band selected galaxies which is significantly younger at 
z ~ 6 than at z ~ 3. 
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